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Abstract. We consider the Abelian model of a complex scalar field coupled to a gauge 
field within the framework of General Relativity and search for cosmological solutions. 
For this purpose we assume a homogeneous, isotropic and uncharged Universe and a 
homogeneous scalar field. This model may be inserted in several contexts in which the 
scalar field might act as inflaton or quintessence, whereas the gauge field might play 
the role of radiation or dark matter, for instance. Particularly, we propose two such 
models: (i) in the first, the inflaton field decays to massive vector bosons that we regard 
as dark- matter; (ii) in the second, due to its coupling to radiation the scalar field is 
displaced from its ground state and drives an accelerated expansion of the Universe, 
playing the role of quintessence. We observe that the equations are quite simplified 
and easier to be solved if we assume a roughly monochromatic radiation spectrum. 
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1. Introduction 

The study of scalar fields in cosmological models has been proven to be very fruitful 
both in inflationary theories and in quintessence models. [TJ El HI 13 E] 

It turns out that the success of these models also raises the need to investigate the 
interaction of these scalar fields with other forms of matter. 

The inflationary theories are based upon the dynamical evolution of a scalar 
field whose coupling to other fields is normally ignored until reheating process takes 
place El IU IZ] • In the first reheating models a phenomenological friction term Tcf) was 
inserted into the equation of motion of the scalar field in order to simulate the energy 
transfer to other forms of matter jH] . The theories of reheating have been improved since 
then, and today they are very more sophisticated, taking under consideration quantum 
effects and many kinds of possible decays of the inflaton to other particles [HJ CDj- 
Anyway, how reheating really occurred is still an object of intense investigation [TTj . 

The very recent measurements of type-IA supernovae ^2] and the analysis of the 
power spectrum of the CMBR provided strong evidence for a present accelerated 
expansion of the Universe ^21 QSj. That offered a new setting for the scalar field 
cosmological models. In these contexts they are named quintessence fields [5], and 
recently many interesting phenomenological models of interacting quintessence have 
been proposed jTo] . 

In the present work a very simple cosmological model is investigated in which a 
complex scalar field interacts with a vector field. We believe that the kind of interaction 
treated in this model is very special because it raises naturally by requiring a local 
£7(1) gauge symmetry. This is the principle beyond the electromagnetic interaction and 
it inspired the development of the other fundamental interactions. According to the 
current view all the forms of interactions are completely determined by certain internal 
symmetry groups [T7j . 

Our starting point will be the following Lagrangian: 

£ ma t = d^(d^)* - v(<m>*) - \f, u f^, (1) 

where 

Dfi = V„ + iqAp, = d^Ay - d v A„. (2) 

The Lagrangian (jJJ is clearly invariant under local U(l) gauge transformations: <f) — > 
exp[iqa(x)](p and A^ — > A^ — d^a{x), and associated with this symmetry is the 
conservation of the four-current 

J^iq[4>{D^Y-^{D^)]. (3) 

Our basic procedure will be to introduce the Lagrangian in the action with 
the Einstein-Hilbert Lagrangian and derive the Einstein's equations for that system 
(Sec. 2). Then we shall investigate solutions to those equations that may be applicable 
in cosmological models. In Sec. 3 we ascribe an interpretation to the contribution of each 
field for the energy density. We also discuss about the flow of energy among the fields. 



Cosmological models of coupled scalar and Abelian gauge fields 



3 



The Fourier expansion of the vector field and spatial average procedures are performed 
in Sec. 4 as a first step towards the application to concrete models. Considerations about 
the possibilities of such applications comes in Sec. 5, and the equations are solved for the 
simplified case of a roughly monochromatic radiation spectrum. Finally, we conclude in 
Sec. 6 with brief remarks and future prospects. All the fields are treated classically. We 
choose units such that c — h—k — 1 and adopt the convention that Greek indices run 
from to 3 whereas Latin indices run from 1 to 3. 



2. Field equations 

As was previously mentioned, this model can describe the interaction of the scalar field 
with the electromagnetic field, provided we identify the gauge vector field with the 
electromagnetic four-potential For the time being, that is what we shall bear in 
mind in our procedure. However, in later sections we will investigate cases in which 
spontaneous symmetry breaking occurs, so it is convenient to keep an open mind in 
the sense that the vector field does not necessarily need to be identified with the 
electromagnetic vector potential. 
Let us consider the action: 



S = J d^X\f^g 



U +D^{D^y-V^(jf)--F IJiV F^ 



(4) 



.16ttG Mrv T ' v " ' 4 
By varying the action with respect to <p*, <fi and A^, we obtain the equations: 

D,D^ = (D.Dyy = F^, v = -J^ (5) 

respectively. By the anti-symmetry of F^ v and the symmetry of R^, eq.([5|U implies 
that J% = -F^w = -R^F»» = 0. 

The Einstein's field equations follow by requiring a stationary action with respect 
to variations of the metric g^ u , and turn out to be: 

Gfiv = -SttG {T^ d + 2 J D (M 0[A, ) 0r - 9^ [D*<f>{D°<f>y - IW)]} , (6) 

where T^ d is the well-known Maxwell energy-momentum tensor of the electromagnetic 
field PI: T- d = F°F av + g^F^/A. 

We are interested in solutions to those equations that describe a homogeneous scalar 
field interacting with a homogeneous and isotropic fluid of incoherent radiation. In that 
purpose, we will stress some points and make some considerations: 

(i) the space-time geometry will be that of a plane Robertson- Walker metric, 
ds 2 = dt 2 - a(t) 2 (dx 2 + dy 2 + dz 2 ); 

(ii) we will assume a homogeneous scalar field, = <f>(t); 

(iii) we will describe an uncharged Universe, such that the charge density J° = 0; 

(iv) the above conditions together with J^-^ = imply that diA 1 = . But that 
does not completely determine the field A^, and we still have the gauge freedom to 
choose A = (Coulomb gauge); 
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(v) the choice A = and the \x = component of eq. ©3 imply that 00* — 4><p* = 0. 
By writing as 0i + 202, we see that the last equality means that 02 — ct'Pii and that 
allows us to define the real field $ such that <3> 2 = 00* = 2 + 2 = (1 + a 2 )0f; 

(vi) by considering the free-field Lagrangian £ f £ ee = (9^0(9^0* — V(00*), we can 
use Noether's theorem to derive the canonical energy-momentum tensor for the scalar 
field, Sft u = 2d( ll (j)d u )(j)* — g^ip^cf — V). Regarding S^ v as the stress-tensor of a 
perfect-fluid, we can make the following identifications: 

p„ = S ° = $ 2 + V(<1> 2 ), p* = —St* = $ 2 - V($ 2 ), (7) 

where p$ and p<j> are the energy density and pressure of the scalar field, respectively; 

(vii) we will assume that T™ d is the energy-momentum tensor of a homogeneous 
and isotropic fluid of incoherent radiation. That means that in a co-moving system 
there is no average net flow of electromagnetic energy in any direction [T9"l |2*U] , so that 
T^ d is diagonal and p r = p r /3, where p T is the radiation pressure and p T is the radiation 
energy density given by p T = T^ 1 . 

In view of the above considerations we can write the system of equations for a(t), 
$(£), and A,i{x^) (bearing in mind that A = 0) as: 

S ^% + P*-g 2 ^$ 2 l, (8) 



H 2 



3 

a 4tiG 



a 3 



[p r + P$ + 3(p r + p$)], (9) 



$ + 3#$ + = q 2 A i A l $, (10) 
2 d§ 



d^d^Ai + HAi = -2q § A { . (11) 

Equations (jH)) and Q refer to the two linearly independent Einstein's field equations 
flOJ), whereas (HHJ) and (|TTJ) follow from ©. 



3. Balance equations and interpretation 

It is well-known from elementary electromagnetic theory ^HJ that the magnetostatic 
energy of a current distribution is: U — | J v J • ArfV^, where J is the current density and 
A is the magnetic vector potential. Hence, it is possible to associate an energy density 
I J • A to a system of currents. Regarding to our system, in which J % = —2q 2 A l Q 2 , the 
current energy density is given by 

PJ = l -AiJ l = -q 2 A % A% 2 . (12) 

Therefore, we can rewrite the Friedmann equation (jHJ) in its familiar form H 2 = 
8nGp/3, and ascribe to it a simple physical interpretation: the total energy density of 
the system, p, is the sum of three forms: 

(i) p r , the energy density of the electromagnetic field, which may also be regarded 
as the radiation energy; 
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(ii) p$, the energy density of the scalar field; 

(iii) pj, the current energy density Moreover, the current pressure is given by 
pj = —pj/3, as can be seen by identifying the energy-momentum tensor with that 
of a perfect fluid: = diag(p, — p, —p, —p), where p = p r +P®+Pj andp = p r +P<s>+P 

Let us now write the balance equations for p T , p$, and pj. First from T^ A . V = 
—F^ U J V it follows the balance equation for the energy density of the radiation field, i.e., 

p I + 3H(p I +p I ) = -J i A i . (13) 

We note that the term — J % Ai = —J ■ E represents the power flux transfered to 
the electromagnetic field due to the currents motion [18 . Next, differentiation of 
p$ = $ 2 + V A ($ 2 ) with respect to time together with (fT0|) lead to 

p # + 3H (p$ + p$) = 2q 2 A i A i <!><5>, (14) 

which is the balance equation for energy density of the scalar field. Finally, from the 
expression for pj, eq. (JEJ), it follows the balance equation for the energy density of the 
currents 

pj + 3H( Pj + pj) = J l M - 2q 2 A t A^<S>. (15) 

Summation of eqs. ((EH), (JHJ) and (fTKj) leads to the balance equation for the total energy 
density of the system 

p + 3H(p + p) = 0, (16) 

naturally not independent of the Friedmann equation (jHJ) and of the acceleration 
equation (JHJ). 

The interactions among the fields might be depicted in the form of a diagram, as 
shown below: 









2q 2 A l Ai<S><l> 




Pr 


< — > 


PJ 


< ¥ 


P$ 



4. Fourier expansion and spatial average procedure 

We now turn to the problem of finding a solution to the system of equations (jH)) through 
|TT|) . Since we will have to deal with a partial differential equation for the field Aj, our 
approach will be to expand it in a Fourier series and make some simplifying assumptions 
in order to make the analysis manageable. All of these assumptions turn out to be quite 
reasonable, and as we shall argue later, are quite well justified. 

The expansion for the vector potential in a triple Fourier series can be written as: 

A i (t,x)=J2Mf(t)e ik -*, (17) 

k 

where the coefficients in the expansion are related by M" k = (Mf)*, since Ai is real. 
From the condition diA 1 = it follows that for each k, kiMf = 0, where hi refers to the 
components of the wavevector k. 

|| It is interesting to note that, due to such barotropic relation (pj = —pj/3), the current energy 
density does not contribute to the acceleration, eq.©, since such a contribution would show up in the 
form pj + 3pj, which is identically zero. 
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The equation satisfied by the coefficients Mf(t) follows from (fTTj). and reads: 



A: 2 



Mf + HMf + — Mf = -2q z $ z Mf. (18) 
a 

Now, we argue that an average procedure is needed in order to describe a 
homogeneous and isotropic radiation field ^HIED]- To that end we define the volumetric 
spatial average [21] of a quantity ip at time t by: 



iff = Tr [ ipy/=gd 3 x, 
V Jv 



(19) 



where V = L 3 is the volume of a fundamental cube. 

Imposing periodic boundary conditions, the summation in expression (JT7jl extends 
over the discrete spectrum of modes of the wavevector k whose components run through 
the values fc, = 27rnj/L, assuming integer values. In the limit that L is large compared 
to the length scale in which homogeneity holds, the passage to the continuum description 
is well justified. 

In the light of this we shall replace the expressions for AiAi and p T = by its 
averaged values in order to make explicit their exclusive dependence on the variable t. 
In this case from 



AiAi 



£M*(t)e"« 



kx 



ik'-x 



1 

2a? 



A^ + —(OiA/lA, - diAjdjAi) 



we get the averaged values 

A-A-{t) = Y,Mf{t)Mf{ty 



(20) 
(21) 

(22) 



2a(t)' 



E 



k L 



(23) 



We remark that the average procedure adopted here, besides making AiAi and p r 
functions of time t only, accounts for the vanishing of the off-diagonal components of the 
energy- momentum tensor T^, which is required for self-consistency since the Einstein 
tensor G^ u is diagonal in our homogeneous and isotropic geometry. 

Now, since Mf{t) is a complex quantity, we will write it in its polar form: 

Mf(t) = P^(t)e iuJ ^\ (24) 

and the equations for the real coefficients P^(t) and ujf(t) can be obtained by inserting 
into (|18|). Hence it follows 



+ Huf + 2^ k 



0. 



1,2,3), 



which can be readily integrated and gives: 
k c} 



UJi = — - 



2 • 



(25) 



(26) 
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where c k is a constant. Moreover 

k 2 



P k + HP* + 



C: 



k2 



2g 2 $ 2 P k 



a* \ ' p k 
In terms of P*(t) (and dropping bars), we have 
A i A i (t) = J2P^t) 2 , 



Pr(t) 



2a(f) ; 



E 

k,i 



a(t) 1 



(27) 



(2? 



(29) 



In short, we have the system of equations (jUJ), (fTU|) . and (J27|) to be solved for 
the functions a(t), $(£) and the infinitely many P*(t) coefficients. Solving for each 
P*(t) means to find out how each frequency band of the radiation spectrum interacts 
and evolves. This task is certainly worth being investigated, and will be performed in a 
future work. Here, for the sake of simplicity we shall solve the problem for the particular 
case in which all the frequencies can be ignored except a specific one, i.e., for the case of 
a monochromatic spectrum. This procedure, to be considered in the next section, will 
make the analysis more tractable and will allow us to gain insight on the qualitative 
dynamical behavior of the model. 



5. Monochromatic spectrum 

For a monochromatic spectrum, the system of equations (jHJ), (fTUj) . and (fTT|) is simplified 
and reads: 

d AttG r , , . 

[p r + p$ + 3Qo r + p<3>)J, (30) 



a 3 

IdV q 2 
2d® ~ ~a 

k 2 



<5> + 3H$ + -— = -^A i A i <&, (31) 



ag(AA) + HdoiAAi) = -Aq^AAi + 4a 2 p r - ^AA, (32) 

where p$ = $ 2 + V, p<s> = $ 2 — V, pj = q 2 AiA^ 2 / a 2 and p r is given by the Friedmann 
equation, namely, 

tt2 8ttG , . 

H = -g-(Pr +P$ +PJ), 

and is related to its pressure by p T = p r /3. 

That system describes a homogeneous and uncharged Universe filled with 
monochromatic radiation coupled to a homogeneous scalar field. Once we have chosen 
an appropriate potential V(Q 2 ) and the initial amount of energy of each field, that 
system will yield a unique solution describing the subsequent evolution of the Universe 
and its constituents. 

Now it is an opportune moment to discuss the choices for the scalar field potential 
V(& 2 ). Up to now we have been identifying the vector field A^ with the electromagnetic 
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Figure 1. The exit of inflation and the beginning of deceleration. 

vector potential. The scalar field potentials V(Q 2 ) that fit in this case are those which 
have a minimum at $ = 0. However, potentials whose lowest state is not at $ = 
might not be suitable as concerns to radiation, since the symmetry will be spontaneously 
broken and the "photons" will become massive. Therefore, we will distinguish between 
two classes of potentials V($ 2 ): the ones with non-degenerate minimum and the rest. 

In what follows we shall investigate the solutions to equations (J3UJ), (j3"Tj) and (|3"2"j) 
in two contexts: reheating (subsection 5.1) and quintessence (subsection 5.2). 

5.1. Reheating and Dark-Matter 

Among the several situations in which this model may be applied, perhaps the most 
straightforward one is the period that succeeds the end of inflation, when the energy 
stored in the inflaton field decays to radiation and reheats the Universe. The most 
suitable potentials V(§ 2 ) in this case are V oc $ 2 or V oc $ 4 . Investigation of eqs. 
(|3U|) to (|32|) in this context does not reveal any new feature than those already known 
from the literature in reheating processes EI]. In fact, this description faces the 
very difficulties that decaying oscillating scalar fields display, namely, that the decay 
products of the scalar field are ultrarelativistic and their energy densities decreases due 
to expansion much faster than the energy of the oscillating field, which goes as cx 1/a 3 . 
Therefore the energy stored in the inflaton field never decays completely, even when the 
coupling q 2 is large. As was remarked in ^U|, reheating can be complete only if the decay 
rate of the scalar field decreases more slowly than Typically this requires either 
spontaneous symmetry breaking or coupling of the inflaton to fermions with < 

Nevertheless, this model also exhibits effects of resonance jUH IHj that may turn 
reheating efficient, although that requires a good amount of fine tuning in the value of k 
(see eq. (J3"2"]) ). that is, the system (|3Dj). (pIT]) and (}3"2"j) would not be the most appropriate 
to allow for resonance behavior, and for a more accurate description the full spectrum 
of radiation must be taken into account. 

Let us now lay radiation aside and allow the photons to gain mass. In that purpose, 
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Figure 2. The growth of the scalar field. 



we will choose the potential V(<& ) = /ie , where p and A are constants. It is clear 
that, as $ rolls down the potential hill towards non-zero values, the photons acquire 
an effective mass. Therefore this choice of potential can describe a scalar field that 
decays to massive vector bosons, which, in turn, can be responsible for dark matter in 
the Universe. To carry that out, we have to abandon our previous interpretation of the 
energy density of the fields which distinguish between pj and p r , and ascribe to the 
massive vector field the energy density pdm = Pr + Pj and pressure pd m = Pr + Pj, where 
the subscript alludes to dark matter. In this redefinition, p r accounts for the kinetic 
part of the energy of the massive photons, while p j owes to the contribution of its mass. 

By numerical analysis we extract the general features of the model. In figures 1 to 
3 we plot the relevant intervals for d(t), and the energy densities of the fields for a 
particular choice of parameters: 



k = 2aoH Q , 

and initial conditions: 

a(0) = a , 

6(0) = o, 

where M p is defined by M p 



q 



0.05# 2 M" 2 , 



A = 10M" 



p = 0.7 p c 



d(0) 

Pr(0) 



V(0) 



k 2 AiA(0)/at, 



0.6p , 



0. 



HqMi 2 . However, the behavior 



'3/(87rG) and p 

displayed is quite general and can be modeled by an adjustment of the parameters and 
initial conditions. It is described as follows: during inflation the scalar field is next to 
$ = 0, where the potential V(Q 2 ) is nearly flat and attains its maximal value. When 
the scalar field begins to depart from $ = and roll down the potential, inflation 
ends and the scalar field decays to massive vector bosons. The decay process is not 
abrupt; rather, there is a smooth transition from a domination of the scalar field to a 
domination of matter. Numerical solutions show that when the decay rate 2q 2 A i A t §& 
of the inflaton (see eq. (|14jl) ceases to be significant, the scalar field attains a roughly 
constant value, $ c , implying an effective mass of the vector bosons md m = V^Q^c (see 
eq. (fTTj)). which, by this time, behave as ordinary matter with p dm ps and pd m scaling 
as oc 1/a 3 , whereas p$ decays faster than oc 1/a 4 . 
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Figure 4. Accelerated expansion driven by the quintessence field. In the small frame: 
end of quintessence era. 



5.2. Quintessence 

By taking another look at the equation of motion for the scalar field: 

... 1 2 

$ + 3H$ + -V = ~AiAi$ 
2 a 1 

we can spot a new term at the right hand side of the equality which is absent in the case 
of a non-interacting scalar field. This term shows up due to the interaction with the 
vector field, and it resembles an elastic force (F = —kx) whose time-dependent "elastic 
constant" is q 2 AiAij a? . This force acts on the motion of the scalar field by drawing 
it towards the position $ = 0. When the scalar field is subject to a potential V(§ 2 ) 
whose minimum lies at $ = the influence of this force will not alter significantly the 
qualitative behavior of the scalar field, but, on the contrary, if $ = is not a minimum, 
there must arise relevant modifications in its motion. 
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Figure 5. The behavior of the field <i>: (i) during the quintessence era; (ii) in the 
subsequent period when it rolls back to its potential minimum and oscillates about it 
(small frame). 

That is what we shall investigate in this section by considering the "Mexican Hat" 
potential 

V{$ 2 ) = A($ 2 - $ 2 ) 2 , where $ 2 = $g 

corresponds to the lowest state. It is clear that the vacuum solution $ 2 = $q and = 
satisfies eqs. (fTUJ) and (fTTf . 

However, in the presence of vector bosons (i.e., if AiAi ^ 0), the ground state for 
the scalar field becomes an unstable position, for if we place $ at $q, it will be pulled 
toward the origin due to the term — q 2 A i A i ^ja l . 

With the following choice of parameters: 

k 2 = 0.7a 2 H 2 , q 2 = 0.5H 2 M~ 2 , A = O.Ol^M" 2 , $ 2 = 0.05M 2 , 

we have solved numerically the system (j30j) . (J3TJ) and with initial conditions 

a(0) = oq, d(0) = a H , $(0) = $ , 

6(0) = 0, Pr (0) = k 2 AiAi(p)/<4, d o (AAi)(0) = 0, 

that corresponds to placing $ in its ground state (p$ = 0) and filling the Universe with 
radiation. In figures 4 to 6 we plot the relevant intervals for a(t), and the energy 
densities of the fields. 

The coupling of the scalar field to radiation removes it from its ground state and 
makes it climb up the potential hill to positions near the origin. That raises its potential 
energy and confers on it a quintessential behavior. During the period that $ remains 
next to $ = the acceleration becomes positive, p r goes as oc 1/a 4 and the vector 
bosons can be, meanwhile, regarded as radiation. This period can be made as long 
as one wishes by a suitable choice of the parameters. But, ultimately, the force that 
keeps $ next to zero becomes too weak and $ rolls back to (one of) its ground state(s), 
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Figure 6. The coupling to radiation transfers energy to the field <I>, which later comes 
to dominate and behaves as quintessence (small frame). Energy densities: scalar field 
- straight line; current field - dotted line; radiation field - dashed line. 



oscillating about it. Consequently, the symmetry is spontaneously broken, the photons 
become massive % and the expansion decelerates again. 

As a final remark we note that, in this model, the positive acceleration of the 
Universe is not due to a small non-zero vacuum energy |22 E]- Here, the vacuum 
energy is indeed zero, but the point is that the Universe is not in the vacuum state. It is 
filled with particles that, due to their interactions, displace the fields from their ground 
states and make them act as dark energy. Motivated by this we suggest that a possible 
explanation for the present acceleration of the Universe might reside in the fact that the 
Universe is (obviously) not in vacuum although its energy density is very small, and the 
interactions among its constituents, which so far have been commonly neglected, might 
be responsible for the dark energy. 

6. Conclusions and Outlook 

To sum up, we remark that the qualitative treatment performed here point out to 
interesting mechanisms of accounting for the dark-matter in the Universe (or at least 
a part of it), and also provides a possible candidate for quintessence. However, the 
quintessence behavior displayed here differs from the usual quintessence models in the 
sense that the scalar field is not tending asymptotically to the minimum of its potential. 
Rather, its minimum is not a stable position and the scalar field is removed from its 
ground state due to its coupling to a background field of vector bosons. + 

All the same, our study is not complete and a quantitative confrontation of the 

% As was pointed out in the last subsection, when the symmetry is broken it is convenient to modify 
the interpretation of the energy densities, and ascribe to the massive photons the energy density 
Pr + q 2 ^AiAi/ 'a 2 , p r being the kinetic contribution to its energy and q 2 Q 2 t AiAi/a 2 the contribution 
due to its mass m — ^/2q&o- 

+ A similar mechanism was employed to produce a short secondary stage of inflation. See |23) and 
references therein. 
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predictions yield by these models to observations is needful. Besides taking into account 
the full spectrum, it is also necessary to address the issue of whether the models treated 
here satisfy the constraints from the anisotropy of the CMBR, the large-scale galaxy 
distribution and the SN-Ia data. That will be performed in a future work. 
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